Pseudogap, non-Fermi-liquid behavior, and particle-hole asymmetry 

in the 2D Hubbard model 
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The effect of doping in the two-dimensional Hubbard model is studied within finite temperature 
exact diagonalization combined with cluster dynamical mean field theory. By employing a mixed 
basis involving cluster sites and bath molecular orbitals for the projection of the lattice Green's 
function onto 2x2 clusters, a considerably more accurate description of the low frequency properties 
of the self-energy is achieved than in a pure site picture. The transition from Fermi-liquid to non- 
Fermi-liquid behavior for decreasing hole doping is studied as a function of Coulomb energy, next- 
nearest neighbor hopping, and temperature. In particular, the self-energy component Ex associated 
with X = (n, 0) is shown to exhibit an onset of non-Fermi-liquid behavior as the hole doping 
decreases below a critical value S c . The imaginary part of Sx(w) then develops a collective mode 
above Ef, which exhibits a distinct dispersion with doping. Accordingly, the real part of Ex(w) 
has a positive slope above Ef, giving rise to an increasing particle-hole asymmetry as the system 
approaches the Mott transition. This behavior is consistent with the removal of spectral weight 
from electron states close to Ef and the opening of a pseudogap which increases with decreasing 
doping. The phase diagram reveals that S c ~ 0.15... 0.20 for various system parameters. For 
electron doping, the collective mode of Ex (w) and the concomitant pseudogap are located below 
the Fermi energy which is consistent the removal of spectral weight from hole states just below Ef- 
The critical doping which marks the onset of non-Fermi-liquid behavior, is systematically smaller 
than for hole doping. 

PACS. 71. 20. Be Transition metals and alloys - 71.27+a Strongly correlated electron systems 



I. INTRODUCTION 

The nature of the metal insulator transition as a func- 
tion of doping is one of the key issues in strongly cor- 
related materials^ Experimental studies of many high- 
T c cuprates reveal a rich phase diagram, with conven- 
tional Fermi-liquid behavior in overdoped metals and 
an anomalous pseudogap phase in underdoped systems 
close to the Mott insulator. One of the most intriguing 
and challenging aspects of the non-Fcrmi-liquid phase 
is the observation of highly non-isotropic behavior in 
momentum spaced Whereas along the nodal direction 
TM well-defined quasiparticles exist, in the vicinity of 
X = (tt, 0) strong deviations from Fermi-liquid behav- 
ior occur. In particular, below a critical doping a pseu- 
dogap appears which becomes more prominent close to 
the Mott insulator. This transition from Fermi-liquid to 
non-Fermi-liquid properties has been widely investigated 
in recent years, and several theoretical models have been 
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Dynamical mean field theoryi^iiLi 8 .^^^ (DMFT) 
provides an elegant and successful framework for the 
description of the correlation induced transition from 
metallic to Mott insulating behavior. 22 The local or 
single-site version of DMFT, however, focusses ex- 
clusively on dynamical correlations which can give 
rise to spectral weight transfer between low and high 
frequencies. To address the momentum dependence 
of the self-energy, it is important to allow for spatial 
fluctuations, at least on a short-range atomic scale. 
For this purpose, several approaches based on cluster 



extensions of DMFT^iS as well as cluster per- 
turbation theory^ 7 - have been proposed. The general 
consensus that has emerged from many studies in this 
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is that scattering processes are indeed much stronger 
close to (tt, 0) and (0, tt) than in other regions of the 
Brillouin zone. Thus, Fermi-liquid behavior first breaks 
down in the antinodal direction and a pseudogap in 
the density of states opens up. In the nodal direction 
between (0, 0) and (tt, tt) Fermi-liquid behavior persists 
and well-defined quasiparticles can be identified. 

In the present work we use exact diagonalization 67 
(ED) in combination with cellular DMFT— (CDMFT) 
to investigate the two-dimensional Hubbard model on 
a square lattice for 2 x 2 clusters. For computational 
reasons, ED has previously been applied to study this 
model at T = i 40 i 44 i 45 i 49 Here, we employ an exten- 
sion to finite temperatures by making use of the Arnoldi 
algorithm 6 — which provides a highly efficient evaluation 
of excited states. Moreover, the cluster ED/DMFT is for- 
mulated in terms of a mixed basis involving cluster sites 
and bath molecular orbitals which allows a very accurate 
projection of the lattice Green's function onto the 2x2 
cluster. 69 Thus, despite the use of only two bath levels 
per cluster orbital (12 levels in total), the spacing be- 
tween excitation energies is very small, so that finite-size 
errors are greatly reduced, even at low temperatures. As 
a result of these refinements, extrapolation from the Mat- 
subara axis yields very accurate self-energies and Green's 
functions at low real frequencies. The same approach 
has recently been used to evaluate the phase diagram 
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of the partially frustrated Hubbard model for triangular 
lattices^ 

The focus of this work is on the transition from Fermi- 
liquid to non-Fermi-liquid behavior for decreasing hole 
and electron doping. In particular, we study how this 
transition varies as a function of Coulomb energy, next- 
nearest neighbor hopping, and temperature. A system- 
atic study of this variation is needed to explore the phase 
diagram of the two-dimensional Hubbard model and has 
to our knowledge not been carried out before. 

The key quantity which exhibits the change from 
Fermi-liquid to non-Fermi-liquid behavior most clearly 
is the self-energy component associated with X = 
(-7T, 0). For hole doping 5 < 15 . . . 20 %, spatial fluctu- 
ations within the cluster give rise to a collective mode 
in the imaginary part of Sx(w) above Ep, in agreement 
with early work for 6 — 0.05 by Jarrell et al^ based 
on quantum Monte Carlo (QMC) calculations within the 
Dynamical Cluster Approximation 26 (DCA). The real 
part of Yix(<jj) then exhibits a positive slope, implying 
removal of spectral weight from electron states close to 
Ep and the opening of a pseudogap in the density of 
states. The evolution of this correlation-induced collec- 
tive mode with decreasing doping leads to a widening 
of the pseudogap until it merges with the Mott gap at 
half-filling. In this region, the density of states acquires 
a very asymmetric shape. At large doping the Fermi 
level is located at a peak in the density of states, while 
for decreasing doping Ep gradually shifts into the pseu- 
dogap, giving rise to a marked particle-hole asymmetry 
in the spectral distributions A(k, uo) due to the reduced 
spectral weight above Ep. Moreover, with decreasing 
doping the pseudogap appears first along the antinodal 
direction before it opens across the entire Fermi surface. 
These results are in excellent correspondence with recent 
angle-resolved photoemission data for Bi2Sr2CaCu20s+5 
by Yang et alQ 

The phase diagram shows that the change from Fermi- 
liquid to non-Fcrmi-liquid behavior is remarkably sta- 
ble, S c w 0.15 . . .0.20, when system parameters, such as 
Coulomb energy, temperature, or second-neighbor hop- 
ping, are varied. For electron doping, the resonance of 
Sjs:(w) is located below the Fermi energy, as expected 
for the removal of hole states just below Ep. The doping 
which defines the onset of non- Fermi- liquid behavior, is 
systematically smaller than for hole doping. Finally, the 
Mott transition induced by electron doping exhibits hys- 
teresis behavior consistent with a first-order transition. 
In the case of hole doping, hysteresis behavior could not 
be identified at the temperatures considered in this work. 
Thus, within the accuracy of our ED/CDMFT approach, 
this transition is either weakly first-order at very low tem- 
peratures or continuous. 

The outline of this paper is as follows. Section II 
presents the main theoretical aspects of our finite T clus- 
ter ED /DMFT approach. Section III provides the results 
2x2 clusters. In particular, we discuss the Mott tran- 
sition, the non-Fermi-liquid properties, the pseudogap, 



electron doping, the phase diagram, and the momentum 
dependence. A summary is presented in Section IV. 



II. CLUSTER ED/DMFT IN MIXED SITE / 
ORBITAL BASIS 

In this Section we outline the finite-temperature ED 
method in the mixed site / molecular orbital basis which 
is employed as highly efficient and accurate impurity 
solver in the cluster DMFT. Let us consider the single- 
band Hubbard model for a two-dimensional square lat- 
tice: 

H = - + H.c.) + UJ2^n n (1) 

(ij)<r i 

where the sum in the first term extends up to second 
neighbors. The band dispersion is given by e(k) = 
— 2t[cos(k x ) + cos(ky)] — 4t' cos(k x )cos(k y ) . In order to 
approximately represent hole-doped cuprate systems, the 
nearest-neighbor hopping integral is defined as t = 0.25 
(band width W = 2). The next-nearest-neighbor integral 
is mainly defined as t' = — Q.3t, but t' = will also be 
considered. The local Coulomb interaction is taken to be 
U = lOt = 2.5 and U = 6t = 1.5. Thus, at half-filling, 
the system is a Mott insulator. (For t' = 0, QMC/DMFT 
calculations for 4-site clusters^ yield U c w 1.4... 1.5, 
in agreement with ED/DMFT results for 2-site and 4- 
site clusters^ These values are consistent with recent 
QMC/DCA calculations for 8-site clusters^ which give 
U c w 1.4... 1.6.) 

Within CDMFT 2 ^ the interacting lattice Green's func- 
tion in the cluster site basis is given by 

Gij(iu> n ) = ^2 i iuJ n + /•* - t(k) - E(iu n )}7. (2) 

k 

where the k sum extends over the reduced Brillouin Zone, 
uj n = (2n + l)7rT are Matsubara frequencies and /i is 
the chemical potential. i(k) denotes the hopping matrix 
for the superlattice and Sy (iw n ) represents the cluster 
self-energy matrix. The lattice constant is taken to be 
a = 1 and site labels refer to 1 = (0,0), 2 = (1,0), 
3 = (0,1), and 4 = (1,1). In this geometry, all diag- 
onal elements of the symmetric matrix Gy are identical 
and there are only two independent off-diagonal elements: 
G12 = G13 = G24 = G34 and G14 = G23. By definition, 
both the lattice Green's function Gy and self-energy Ey 
have continuous spectral distributions at real lo. Only 
the paramagnetic phase will be considered here. 

It is useful to transform the site basis into a molecu- 
lar orbital basis in which the Green's function and self- 
energy become diagonal. The orbitals are defined as: 
<t> x = (|1) + 12) + 13) + 14»/2, 4>2 = (|1> - |2) |3> + |4»/2, 
03 = (|1) + |2) - |3) - |4»/2, 04 = (|1) - \2) + |3) - |4))/2. 
We refer to these orbitals as T, M and X, respectively, 
where X is doubly degenerate. The Green's function el- 
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FIG. 1: (Color online) Total density of states p(uj) and molec- 
ular orbital components p m (w) for four-site clusters of square 
lattice. For clarity, the molecular orbital components are di- 
vided by n c = 4. ui — defines the Fermi energy for half- 
filling. At 14 % hole doping the van Hove singularity is shifted 
from uj = -0.13 to -0.09. 



cmcnts in this basis will be denoted as G m (iu) n ), where 

Gr = G\ = Gn + 2Gi2 + G14 
Gm = G2 = Gn — 2G12 + G14 (3) 

X — — Lt4 — Lrll — <J"14- 

An analogous notation is used for the self-energy. Similar 
diagonal representations of G and £ have been used in 
several previous worW 9 i 40 i 53 i 54 i 58 i 59 i 69 i 72 

Figure [T] illustrates the uncorrelated density of states 
components in the molecular orbital basis, where 
p m (w) = — ilmG m (oj) for S = 0, and we denote 
Pr = Pi, Pm = P2, Px = P3A- The average or local 
density is p av = (pr + Pm + 2px )/4. Note that all molec- 
ular orbital densities extend across the entire band width. 
Nevertheless, only px contains the van Hove singularity, 
while pr and pu are roughly representative of the spec- 
tral weight near k = (0,0) and k = (n, n), respectively. 
Hole doping shifts the van Hove singularity towards Ep, 
whereas electron doping moves this singularity away from 
Ep . 

A central feature of DMFT is that, to avoid double- 
counting of Coulomb interactions in the quantum impu- 
rity calculation, the self-energy must be removed from 
the small cluster in which correlations are treated explic- 
itly. This removal yields the Green's function 



Go(iuj n ) = [G(iuj n ) 1 + T.(iLo n )]' 



(4) 



which is also diagonal in the molecular orbital basis. 

For the purpose of perfoming the ED calculation we 
now project the diagonal components of Ga(iuj n ) onto 
those of a larger cluster consisting of n c — 4 impurity 
levels and nb = 8 bath levels. The total number of levels 
is n s — n c + rib — 12. Thus, 



Go :m (iijJ n ) 
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FIG. 2: (Color online) Projection of lattice Green's func- 
tion components Go :m (iui n ) onto cluster consisting of four 
impurity levels and eight bath levels, for U = 2.5, p — 0.5, 
T = 0.01; t' — — 0.3t. Upper panel: ImGo.m, lower panel: 
ReGo,m- Continuous curves: diagonal elements of lattice 
Green's function, Eq. Q; circles: approximate expression, 
right-hand side of Eq. (f5} • 



where e m denote the molecular orbital levels, efc the bath 
levels, and V m k the hybridization matrix elements. The 
incorporation of the impurity level e m ensures a much 
better fit of Go )m (*^n) than by projecting only onto bath 
orbitals. 

Assuming independent baths for the cluster orbitals, 
each component Go tm (iu> n ) is fitted using five parame- 
ters: one impurity level e m , two bath levels and two 
hopping integrals V m k- For instance, orbital 1 couples 
to bath levels 5 and 9, orbital 2 to bath levels 6 and 10, 
etc. For the three independent cluster Green's functions, 
we therefore use a total of 15 fit parameters to represent 
Go(iuj n ). This procedure provides a considerably more 
flexible projection than within a pure site basis. Since 
for symmetry reasons all sites are equivalent one would 
have in this case only four parameters (without includ- 
ing a level at the cluster sites). Thus, the molecular or- 
bital basis allows for 11 additional cross hybridization 
terms as well as internal cluster couplings (see below). 
In addition, it is much more reliable to fit the three inde- 
pendent molecular orbital components Go^ m (iuj n ) than a 
non-diagonal site matrix Go^j(itu n ) with only 4 parame- 
ters. 

Figure [5] illustrates the projection of the lattice Green's 
function Go(iui n ) onto the cluster for U = 2.5 and 
p: = 0.5, which corresponds to about r5 = 0.08 hole dop- 
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ing. Projections of similar quality are achieved at other 
Coulomb energies and chemical potentials. 

The evaluation of the finite temperature interacting 
cluster Green's function could in principle also be car- 
ried out in the molecular orbital basis. The Coulomb 
interaction then becomes a matrix containing many inter- 
orbital components. This step can be circumvented by 
using a mixed basis consisting of cluster sites i and bath 
molecular orbitals k. Thus, the diagonal 8x8 subblock 
hb = {tk$kk') representing the bath levels remains un- 
changed, but the diagonal 4x4 cluster molecular orbital 
submatrix now becomes nondiagonal in the cluster site 
basis. The transformation between sites i and orbitals m 
is given by 



(G) 



In this mixed basis, the site subblock of the cluster Hamil- 
tonian becomes 




he = 




(7) 



with e = (ei + e 2 + 2e 3 )/4, r = (d + e 2 - 2e 3 )/4, and 
t' = (ei — ea)/2. Note that the hopping elements t and 
t' of the original lattice Hamiltonian do not appear since 
they are effectively absorbed into r and t' via the molec- 
ular orbital cluster levels e m which are adjusted to fit 
Go } m(iw n ) ■ Evidently, the procedure outlined above not 
only includes hopping between cluster and bath. It also 
introduces three new parameters within the 2x2 clus- 
ter: e, r, and r'. In the mixed basis, the hybridization 
matrix elements V m k between cluster and bath molecu- 
lar orbitals introduced in Eq. |(5J) are transformed to new 
hybridization matrix elements between cluster sites i and 
bath orbitals k. They are given by 



V( k = (TV) ik = 



(8) 



Thus, the upper right 4x8 submatrix containing the clus- 
ter / bath hybridization matrix elements is transformed 
from 



(9) 
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(10) 



The single-particle part of the cluster Hamiltonian then 
reads 



h 



h c V 
V n h h 



(11) 



Adding the onsite Coulomb interactions to this hamil- 
tonian, the non-diagonal interacting cluster Green's func- 
tion at finite T can be derived from the expressio n 73 ' 74 



E,, — E.. 



H4xIm)(m|c jct H 

En - E v + iuj n 



(12) 



where E v and \i>) denote the eigenvalues and eigenvectors 
of the Hamiltonian, (3 = l/T and Z — ^2, v exp(— j3E u ) 
is the partition function. At low temperatures only a 
small number of excited states in a few spin sectors con- 
tributes to Gfy They can be efficiently evaluated using 
the Arnoldi algorithm. 6 \ The excited state Green's func- 
tions are computed using the Lanczos procedure. Further 
details are provided in Rcf. 73 . The non-diagonal elements 
of Gfj are derived by first evaluating the diagonal com- 
ponents Gfo and then using the relation 



Since G c \ = G C L this yields: 



G, c + G 



j i 



G 



j j • 



G'i 



^ ii ii 1 ' 



(13) 



(14) 



The interacting cluster Green's function Gfj satisfies the 
same symmetry properties as Gy and Go,,j. It may 
therefore also be diagonalized, yielding cluster molecu- 
lar orbital components G^. The cluster molecular or- 
bital self-energies can then be defined by an expression 
analogous to Eqs. ((3]): 



(15) 



At real u>, these cluster self-energy components, just like 
Gq m and G^, have discrete spectral distributions. 

The key assumption in DMFT is now that the impurity 
cluster self-energy is a physically reasonable representa- 
tion of the lattice self-energy. Thus, 



(16) 



where, at real frequencies, S m is continuous. 

In the next iteration step, these diagonal self-energy 
components are used as input in the lattice Green's func- 
tion Eq. @ , which in the molecular orbital basis is given 
by 



G m (iu n ) = ^2 [i^n + /i - Tt(k)T 1 - T,(iu> n )] 



-l 

rnrn 



(17) 
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where T is the transformation defined in Eq. (J6j) . Thus, 
except for the diagonalization which is carried out in the 
mixed site / molecular orbital basis, all other steps of 
the calculational procedure are performed in the diago- 
nal orbital basis. Note that Ti(k)T _1 is not diagonal 
at general k points. As a result, all orbital components 
of £(ia>„) contribute to each G m (ico n ). This feature of 
CDMFT differs from DCA where one has a one-to-one 
relation between H m (iuj n ) and G m (iu>n)'^ 

G» CA (iLj n ) = k + M - e(k) - E^CtWn)]- 1 (18) 

where k m labels the m th patch of the Brillouin zone. 

The largest spin sector for n s = 12 is »f = = 6 with 
dimension AT = 853776. The interacting cluster Hamil- 
tonian matrix h is extremely sparse, so that only about 
20 non-zero matrix elements per row need to be stored. 
Since the Arnoldi algorithm requires only operations of 
the type hu — v, where u, v are vectors of dimension N, 
the procedure outlined above can easily be parallelized. 
At temperatures of the order of T = 0.005 . . . 0.02, one it- 
eration takes about 15 to 60 min on 8 processors. Except 
near the Mott transition, 5 to 10 iterations are usually 
required to achieve self-consistency. 

We conclude this section by pointing out that, once 
iteration to self-consistency has been carried out, a 
periodic lattice Green's function may be constructed 
from the cluster components in Eq. ^ by using the 
superposition^ 

1 4 

ij=l 

= Gn(icj n ) + Gi 2 (iLu n )[cos(k x ) + cos(fcy)] 
+ Gi4(iLu n )cos(k x )cos(k y ) (19) 

At high-symmetry points, this definition coincides with 
the diagonal elements introduced in Eq. (j4|). Thus, 
Gr(iun) = G((0,0),iu n ), G M {iu n ) = C7((7r,7r),iw n ), 
and G x (iuj n ) = G((tt, 0), iu) n ) = G((0, tt), iu) n ). At 
k = (it/2, 7t/2), G coincides with the onsite Green's func- 
tion Gii = (G r + G M + 2G x )/4. 

III. RESULTS AND DISCUSSION 

A. Mott Transition 

Figure [3] shows the occupancies of the cluster molec- 
ular orbitals T, M and X as functions of chemical po- 
tential. The average occupancy per site (both spins) is 
n = (nr + um + 2nx)/2 = 1 — 5, where 5 is the hole 
doping. As revealed by the spectral distributions dis- 
cussed below, the Mott transition occurs at rj 0.7, 
where the X orbital becomes half-filled, whereas nr and 
tim approach 0.25 and 0.75, respectively. Thus, all three 
orbitals take part in the transition. This result is con- 
sistent with previous ED/DMFT calculations^ for 2-site 
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FIG. 3: (Color online) Upper panel: Occupancies of cluster 
molecular orbitals (per spin) as functions of chemical potential 
IM, for U = 2.5, T = 0.01, t' = -0.3t. The Mott transition 
occurs at about fi m 0.7, indicated by the vertical bar, where 
nx — * 0.5, nr — * 0.75, and hm — * 0.25. Lower panel: Average 
double occupancy per site as a function of /i. 

and 4-site clusters in the limit t' = 0, and with recent 
QMC results^ for a minimal 2-site cluster DCA version, 
where hole doping takes place at about the same rate for 
both inner and outer regions of the Brillouin zone. These 
trends differ, however, from results for an 8-site contin- 
uous time (CT) QMC/DCA calculation^ which reveals 
initial doping primarily along the nodal direction, while 
near X the occupancy for small S remains at the same 
value as in the Mott insulator. Evidently, 2-site and 4- 
site cluster DMFT approaches do not provide sufficient 
momentum resolution to allow for k-dependent doping. 

The lower panel of Fig. 3 shows the average double 
occupancy per site. We have calculated these occupan- 
cies both for increasing and decreasing chemical potential 
without encountering hysteresis behavior for T > 0.005. 
The Mott transition induced by hole doping is therefore 
weakly first-order at even lower temperatures, or contin- 
uous. This result differs from the case of electron doping 
discussed farther below, where n m ([x) as well as d occ (ii) 
readily show hysteresis. 

To illustrate the Mott transition in the limit of half- 
filling, we show in Fig. [¥] the spectral distributions ob- 
tained from the interacting cluster Green's function: 
A m (uj) = -(1/tt) ImG^(w + ij), where 7 = 0.02. These 
spectra can be evaluated without requiring analytic con- 
tinuation from Matsubara to real frequencies. The total 
density of states per spin is given by A{ui) = [Aj*(w) + 
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FIG. 4: (Color online) Cluster spectral distributions for var- 
ious chemical potentials corresponding to hole doping S — 
0.12, 8 = 0.08, and 8 = 0; U = 2.5, T = 0.01; broadening 
7 = 0.02. Solid red curves: 4xH; short-dashed blue curves: 
Ar(u}), long-dashed green curves: Am(ui). 



FIG. 5: (Color online) Doping induced spectral weight W+(S) 
above Ef = up to about ui — 0.9, i.e., in the lower part 
of the main gap. A small constant weight is subtracted to 
account for the artificial broadening of the spectral peaks. 
This constant is chosen so that W+ = at /J, = 0.7, 5 = 0. The 
short-dashed blue line defines 8 and the long-dashed green line 
8 + rid , where rid is the average double occupancy shown in 
Fig. 3. 



Our ED/DMFT cluster calculations arc in excellent 
agreement with this picture, as illustrated in Fig. [5] which 
shows the integrated spectral weight per spin induced 
just above Ep, This weight is denoted here as W+(5). 
The initial slope of W+ is seen to be well represented 
by 5, confirming the scenario discussed above. At finite 
doping W+(5) becomes even larger than S + rid, where rid 
is the double-occupancy shown in Fig. [3] These results 
differ from those for t' = and T — obtained by Sakai 
et al£^, who found W + (5) w 8 + rid up to about 14 % 
hole doping. 

Upon closer inspection, the spectral distributions 
shown in Fig. [4] at finite doping reveal a pseudogap close 
to Ep which will be discussed in more detail in the follow- 
ing subsections. As shown below, this pseudogap is inti- 
mately related to the non-Fermi-liquid properties which 
are evident in the X component of the self-energy. 



Am(oj) + 2Ax(w))/4:. All cluster molecular orbitals con- 
tribute to the spectral weight near the Fermi level in the 
metallic phase for 8 > 0, and to the upper and lower 
Hubbard bands in the Mott phase at 8 = 0. 

The evolution of these spectra as a function of dop- 
ing supports the picture conjectured long ago by Eskes 
et alS^ Upon hole doping, spectral weight is transfered 
from the upper and lower Hubbard bands to states just 
above Ep, in the lower part of the Mott gap. Since the 
spectral weight (per spin) of both Hubbard bands ini- 
tially decreases like (I — S)/2, the states induced just 
above Ep have weight 8 (see also RefJ£). This scenario 
is a remarkable consequence of strong dynamical correla- 
tions and differs fundamentally from the one in ordinary 
semiconductors, where states induced in the gap have 
weight 8/2 per spin for total doping 8. 



B. Non-Fermi-Liquid Properties 

We now discuss the low- frequency variation of the clus- 
ter self-energy which is strikingly different for the differ- 
ent cluster molecular orbitals. Fig.[S]shows the imaginary 
parts of E m (zw„), Eq. (|I6p . for chemical potentials ri in 
the range from to 24 % hole doping. The T orbital, 
approximately representative of the center of the Bril- 
louin zone, exhibits the weakest self-energy, ft is nearly 
independent of doping and Fermi-liquid-like, with only a 
moderate effective mass enhancement. Em changes from 
Fermi-liquid behavior at large doping to nearly insulat- 
ing behavior ~ l/iu> n close to the Mott transition. At 
small finite doping, it reveals strong effective mass en- 
hancement. Both fm Er and fm Em extrapolate to very 
small finite values in the limit w„ — > 0, except near the 
Mott transition, fn striking contrast to these orbitals, 
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FIG. 6: (Color online) Imaginary part of self-energy molec- 
ular orbital components E m (ioj n ) as functions of Mat- 
subara frequency for various chemical potentials: /j, — 
0.2, 0.3, 0.4, 0.5, 0.55, 0.6, 0.65 (from top to bottom) cor- 
responding to 5 = 0.24, 0.18, 0.12, 0.08, 0.05, 0.03, 0.01, 
respectively; U = 2.5, T = 0.01. 



ImEx(iu„) exhibits a finite onset in the low- frequency 
limit once the doping is smaller than about 20 % (see 
expanded scale in Fig. [7]) . The onset is largest at about 
fi = 0.55, corresponding to S = 5 %. At smaller dop- 
ing (larger //), i.e., very close to the Mott transition, it 
diminishes again. 

In addition to the low- frequency onset of ImT, x(iuj n ) 
which gives rise to reduced quasiparticle lifetime, the 
non-Fermi-liquid behavior also leads to a characteristic 
flattening of ImS^(iw n ), which induces a sharp reso- 
nance in ImEx(u) at small positive frequencies. As will 
be discussed in the next subsection, it is this resonance 
that is responsible for the pseudogap in the density of 
states. 

Similar results are obtained at lower temperature, T = 
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FIG. 7: (Color online) Im Ex(ia)n) for various hole doping 
concentrations on expanded scale for T — 0.01 (upper panel) 
and T = 0.005 (lower panel). The long-dashed green curves 
near 5 « 18 % denote the approximate onset of non-Fermi- 
liquid behavior which grows until 8 decreases to about 5 %. 



0.005, as shown in the lower panel of Fig. [71 Again, the 
largest deviation from Fermi-liquid behavior is found for 
ImEx(iw„) at about 6 % doping. The onset of non- 
Fermi-liquid properties occurs at slightly smaller doping 
than for T = 0.01. The results shown in Figs. [5] and [7] 
are consistent with the T = ED/CDMFT calculations 
by Civelli et aZiP- 

To illustrate this non-Fermi-liquid behavior of T*x in 
more detail, we compare in the upper panel of Fig. [5] the 
low- frequency limits "f m = — ImE m (iw„ — > 0) as func- 
tions of chemical potential. These values were found to 
be nearly the same for a linear extrapolation from the 
first two Matsubara points and for a quadratic fit us- 
ing the first three points. At fi > 0.3 or 5 < 0.18, -fx 
increases strongly, indicating the onset of a non-Fermi- 
liquid phase. The lower panel shows the variation of -fx 
with doping for T = 0.01 and T = 0.005. At lower T, the 
onset of non-Fermi-liquid behavior is seen to be slightly 
sharper and to shift to slightly lower S. 

At finite temperature, a sharp transition between 
Fermi-liquid and non-Fermi-liquid phases is not to be 
expected. According to the detailed temperature vari- 
ation of the self-energy of the two-dimensional Hubbard 
model studied recently by Vidhyadhiraja et al. e3 within 
QMC/DCA for 4 x 4 clusters (U = 1.5, = 0), a quan- 
tum critical point exhibiting marginal Fermi-liquid be- 
havior was found at 5 C w 15 % doping, with Fermi- 
liquid behavior at larger 6 and a pseudogap phase at 
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FIG. 8: (Color online) Upper panel: Low-frequency damping 
rates 7 m = — ImE m (itj n — > 0) as functions of chemical poten- 
tial. Solid red curves: X orbital, short-dashed blue curves: 
F orbital, long-dashed green curves: M orbital. The dot- 
ted curve denotes the doping 8 (same scale as 7 m ); U = 2.5, 
T = 0.01. For n > 0.2 ... 0.3 or 8 < 0.18 . . . 0.20, jx increases 
strongly, indicating the onset of non-Fermi-liquid behavior. 
Lower panel: comparison of 7x as a function of doping for 
T = 0.01 (solid red circles) and T — 0.005 (empty blue cir- 
cles). The vertical bar denotes the approximate location of 
the transition from Fermi-liquid to non- Fermi-liquid behavior. 



6 < S c . At T = 0.01, the T J 5 phase diagram indicates 
a crossover region of about S = 0.15 ± 0.02 between these 
phases. Assuming a crossover region of similar width, 
i.e., AS ~ 0.04, the results shown in Fig. [5] are consistent 
with those of Ref. 63 . Thus, for U = 2.5, = -0.3t, the 
Fermi-liquid and non- Fermi-liquid phases seem to be sep- 
arated by a quantum critical point at S c rs 18 ... 20 %, 
with marginal Fermi-liquid behavior for T > 0. 

Figure [5] compares the low- frequency damping rate of 
the X orbital as a function of doping for U — 2.5 and U = 
1.5. The upper panel shows the results for t' = — 0.3t, the 
lower panel for t' = 0. The case U = 1.5, t' = suggests 
a critical doping S c » 0.15 ± 0.02, in agreement with the 
results of RefJ£. As is to be expected, at smaller U S c is 
smaller than at large U, since the Fermi liquid properties 
are stabilized. A similar trend occurs as t' is shifted from 
t' = — Q.3t to t' = 0. Nevertheless, despite the large 
variations in U and t', the critical doping separating the 
Fermi-liquid and non-Fermi-liquid phases is remarkably 
stable and occurs in the range S c w 0.15... 0.20, i.e., 
close to the optimal doping concentrations found in many 
high-T c cuprates. 



FIG. 9: (Color online) Low- frequency damping rate of X or- 
bital, 7x = — Im'Ex (it±>n — » 0) as a function of doping for 
various Coumb energies and next-nearest-neighbor hopping 
energies; T = 0.01. Upper panel: t' = — Q.St, lower panel: 
t' — 0. Solid red circles: U = 2.5, empty blue circles: U = 1.5. 



For a more detailed comparison with the results of 
Refi^ 3 -, we show in Fig. [Til] the variation of ImSx(»w„) 
with chemical potential for U = 1.5, t' = and T — 
0.01. These values of fj, correspond to dopings in the 
range S = 0.27... 0.05. Although the overall mag- 
nitude of Im Ex is much smaller than in Fig. [7] for 
U = 2.5, t' = — 0.3t, there is again a clear separa- 
tion between doping larger than S c rj 0.17 exhibiting 
Fermi-liquid behavior, and smaller doping with char- 
acteristic non- Fermi-liquid features. The lower panel 
shows the comparison of the approximate quasiparticle 
weight, Zx = 1/[1 — ImSx(!Wo)/wo], derived from the 
ED/DMFT results in the upper panel, with the corre- 
sponding QMC/DCA values taken from Fig. 1 of Ref. 63 . 
For S > 0.15 the agreement is very good. (Note that for 
Zx > 0.5, — ImSx(iwo) is less than lu — 0.031.) At 
smaller doping, the difference becomes larger, presum- 
ably because of the finer momentum resolution achieved 
for the 4x4 cluster in Refi^. 

To analyze the difference between CDMFT and DCA 
for 2x2 clusters, with identical system parameters, we 
compare in Fig. [TJJ the low-frequency damping rate of 
Yix(ito n ) as a function of doping. Evidently, the different 
relations between self-energy components E m and lattice 
Green's function G m in these two schemes give rise to 
changes on a quantitative level. Nevertheless, both ap- 
proaches predict a transition from a Fermi-liquid phase 
at hole doping larger than about 20 % to a non-Fermi- 
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FIG. 10: (Color online) Upper panel: imaginary part 
of Ex(i^n) as a function Matsubara frequency for fi — 
0.2... 0.55 in steps of 0.05 (from top), for U = 1.5, t' = 
0, T — 0.01. fi = 0.3 (long-dashed green curve) corre- 
sponds to 8 = 0.18 and approximately marks the transition 
from Fermi-liquid to non-Fermi-liquid behavior. Lower panel: 
Solid red circles: Zx = 1/(1 — ImEx(*<x>o)/wo) derived from 
ED/CDMFT results in upper panel; empty blue circles: anal- 
ogous QMC/DCA results from Fig. 1 of Ref. 63 for T = 0.014. 



liquid phase at small doping. Surprisingly, the transition 
is less sharp in DC A than in CDMFT. The reason for 
this difference might be that, in contrast to CDMFT, the 
momentum patches of the Brillouin zone are not coupled 
in the evaluation of the DCA lattice Green's function 
(see Eq. (|18p). It might therefore be necessary in DCA 
to treat larger clusters (such as 8 sites^ or 16 sites^) 
in order to obtain a sharper Fermi-liquid to non-Fermi- 
liquid transition. A slower convergence with cluster size 
in DCA is also found for the critical Coulomb energy at 
half filling^^ 



C. Pseudogap 

The non- Fermi- liquid properties of Yix(iu>o) manifest 
themselves not only in the enhanced low-frequency damp- 
ing rate discussed above, but also in the flattening of 
ImEx(iWn) which can be identified as the origin of the 
pseudogap in the density of states. Narrow gaps near Ep 
below the critical doping are already evident in the clus- 
ter spectra shown in Fig. [U Fig. [T2] shows these spectra 
on an expanded scale for S = S c « 0.18 and 5 = 0.03. 
While near critical doping the density of states is Fermi- 
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FIG. 11: (Color online) Comparison of low-frequency damp- 
ing rate of X orbital, "fx = — Im Ex (iuo n — » 0) as a function of 
doping for CMDFT (solid red circles) and DCA (empty blue 
circles); U = 2.5, T = 0.01. 



liquid- like, with a sharp peak at Ep, smaller hole doping 
leads to a very asymmetric density of states, with a pseu- 
dogap of magnitude A w 4t 2 /U = 0.1 right above Ep. 
The molecular orbital analysis of these spectra reveals 
that this pseudogap is associated entirely with the Ax (w) 
contribution, i.e., with scattering processes involving mo- 
menta close to (ir, 0) and (0, ir). With decreasing doping, 
the peak at Ep seen for 5 w S c shifts downwards, so 
that the Fermi level gradually moves into the pseudo- 
gap. At the same time, the pseudogap becomes wider 
and the spectral weight above Ep is reduced until the 
transition to the Mott phase occurs at half filling. (The 
peak at u> « 0.25 for 8 = 0.18 is due to the discreteness of 
the cluster spectra and is not related to the pseudogap. 
The actual pseudogap at this large doping is vanishingly 
small; see analysis of self-energy below). 

Note that the peak at Ep for S — S c is also compatible 
with marginal Fermi-liquid behavior. Finite-size effects, 
however, do not permit a clear distinction between Fermi- 
liquid properties below the first Matsubara frequency lvq 
and genuine marginal Fermi-liquid behavior at 6 C J& 

The lower panel of Fig. [T^] shows the corresponding 
spectra derived from the lattice Green's function com- 
ponents G m (iuj n ), Eq. (|17p . via extrapolation to real lu. 
Thus, A{uj) = -±Im[G r (w) + G M (u) + 2G x {u)]/A. We 
use here the routine ratintS^ Nearly identical spectra are 
obtained via Pade extrapolation. About 400 . . . 600 Mat- 
subara points are taken into account for the energy win- 
dow — 1 < u) < 1, and the same broadening is assumed 
(7 = 0.02) as in the cluster spectra shown in the up- 
per panel. As a result of the accurate self-energies and 
Green's functions along the Matsubara axis, the extrapo- 
lation to low real u> is highly reliable. The lattice spectra 
confirm the trend observed in the cluster spectra: At 
fi = 0.3, S = 0.18, the density of states has a peak very 
close to the Fermi level, while for [i = 0.6, S = 0.03, Ep 
lies in a pseudogap of about the same width as in the 
cluster data. The lattice spectra A(oj) can also be cal- 
culated by first extrapolating the self-energy components 
Yi m (iuj n ) to real frequencies and then using Eq. (|17p at 
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FIG. 12: (Color online) Upper panel: Cluster spectral distri- 
butions A(us) = |[A r (w) + A m {lu) + 2A x (co)] for U = 2.5, 
t' = -0.3t, T = 0.01. Solid red curve: = 0.3 (hole dop- 
ing 8 — 0.18) with peak at Ef = 0; dashed blue curve: 
(i = 0.6 (S — 0.03) exhibiting a pseudogap above Ef; lower 
panel: analogous lattice spectra obtained via extrapolation of 
Green's function components, Eq. (jXTJ) . 



real us. The results are fully consistent with the spectra 
derived via extrapolation of G m (ius n ). 

The pseudogap seen in Fig. [T2] for 5 = 0.03 is reminis- 
cent of the pseudogap obtained in the two-band model 
within local DMFT above the first Mott transition. 77 
Once the electrons in the narrow subband are Mott lo- 
calized, an effective two-fluid system is realized in which 
the Coulomb interaction with the remaining conduction 
electrons generates deviations from Fermi-liquid behav- 
ior, in particular, the finite lifetime associated with the 
low-frequency limit of Im£(iw n ), and the characteristic 
flattening of this function which gives rise to a pseudo- 
gap at real This two-band model exhibits a quantum 
critical point when the pseudogap turns into the Mott 
gapiS. It would be interesting to inquire whether the 
present cluster picture of the single band model could 
be mapped onto this two-band model. The spatial de- 
grees of freedom in the cluster would then play the role 
of the inter-orbital fluctuations in the two-band model. 
Since at small hole doping a sizable number of electrons 
is Mott localized their spins act as scattering centers for 
the remaining electrons whose self-energy then exhibits 
deviations from Fermi-liquid behavior. 

To illustrate the effect of non-Fermi-liquid behavior 
on the self-energy at real us, we show in Fig. [T3] the 
low-frequency variation of ImSj(w) obtained from the 
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FIG. 13: (Color online) Self-energy £x(w) obtained via ex- 
trapolation to real frequencies for several hole doping concen- 
trations (broadening 7 = 0.005). Upper panels: — ImEx and 
Re Ex ; U = 2.5, T = 0.005. The outer intersections of Re Ex 
with the straight lines us + n — tk yield the approximate width 
of the pseudogap A. Lower panel: Pseudogap A as a func- 
tion of hole doping derived from self-energies Ex (w) shown 
in upper panels. 



lower panel of Fig. [7] via extrapolation to real us. Typi- 
cally, at these low frequencies we use the first 100 . . . 400 
Matsubara points and evaluate Im Sx at o> + £7, with 
7 = 0.005. Although the details of the resulting spectra 
differ slightly, the important qualitative features near Ep 
are very stable. Spectra derived via Pade extrapolation 
are very similar. 

As can be seen in Fig. [T31 at large hole doping 
-ImEx(u) has a minimum at Ep and varies quadrati- 
cally at small us, as expected for a Fermi liquid. Damping 
in this range is very weak. Nevertheless, even for S > 0.17 
a small peak in -ImEx(w) is found at about lo = 0.05 
above Ep, indicating that electrons added to the system 
just above Ep have a reduced lifetime. With decreasing 
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doping, this feature grows into a strong resonance which 
eventually dominates the low-frequency properties. The 
minimum of -ImEj((w) is then shifted slightly below 
Ep and a second minimum appears above Ep. Moreover, 
this resonance shows a dispersion as a function of doping. 
It first shifts downwards from u> « 0.05 to 0.02 and then 
disperses again upwards to about u> ps 0.05. ReEjf (u>) is 
seen to exhibit a positive slope at the resonance which is 
consistent with Kramers-Kronig relations. This implies 
that spectral weight is removed from the resonance region 
where correlation induced damping is large. 

The evolution of the resonance in ImEj(w) with dop- 
ing is one of the main results of this work and has to our 
knowledge not been discussed previously. A weak reso- 
nance in Im at 5 % doping was also found by Jarrcll 
et alM within QMC/DCA for n c = 4. The fact that this 
resonance is much stronger in the present results might be 
related to the faster convergence of CDMFT with cluster- 
size (see discussion of Fig. [TT|) . A resonance in Im E(u;) is 
also obtained in the spectral weight transfer model pro- 
posed by Philipps et alJ£- In this scheme, however, the 
resonance is located at lo — independently of doping. 

The outer intersections of ReY<x(<+>) with lo + fi — 6k 
provide the approximate width of the pseudogap A in 
the spectral distribution. The central intersection does 
not yield any peak because of the short lifetime in this 
frequency range. The new minima of -ImEx(w) below 
and above the resonance are consistent with the spectral 
peaks just below and above Ep, as seen in the results 
for (5 = 0.03 in Fig. [12] For increasing hole doping, the 
resonance of ImSx(w) becomes weaker so that for S > 
0.17 there are no longer three intersections of u) + /j, — 
€k with ReSjs:(cj). The pseudogap then vanishes. At 
smaller doping, the peak in -ImEj(u) grows and the 
pseudogap gets wider. This trend, however, is superceded 
by the reduction of spectral weight above Ep as the Mott 
transition at half-filling is approached. 

The lower panel of Fig. [T3] shows the approximate 
width of the pseudogap A. We use here the outer inter- 
sections of Re Y>x (w) with the lines lo + fi — to define 
the magnitude of A, where tu is chosen so that lo + /i — 
passes through the inflection point in the region of the 
maximal positive slope of ReSx(^). Other values of 
yield similar values of A. In the spectral distributions, 
this definition of the pseudogap roughly corresponds to 
the peak-to-peak separation of spectral weight near the 
gap. Systematically smaller values of A are obtained, for 
instance, if the width of the gap half-way between the 
minimum of A(lo) and the neighboring maxima is chosen 
as definition. At 5 > 0.17, the definition used above no 
longer yields a pseudogap and the system turns into an 
ordinary Fermi liquid. 

The doping dependent resonance in ImEx(^) and the 
concomitant opening of the pseudogap are consistent 
with recent angle-resolved photoemission (ARPES) data 
by Yang et alZ^ According to the upper panel of Fig. [TBI 
for S > 0.17 the quasiparticle damping is symmetric for 
electron and hole states. Below this doping, the lifetime 



of electron states above Ep is much shorter than that 
of hole states below Ep, giving rise to the opening of 
the pseudogap above Ep and the striking particle-hole 
asymmetry observed in the data. Moreover, the results 
shown in Fig. [13] are specific to the (tt, 0) component 
of the self-energy and are absent in Sr(^) and Sm(u). 
Thus, the particle-hole asymmetry and pseudogap above 
Ep are momentum dependent features which are most 
pronounced in the antinodal region, but weak or absent 
along the nodal TM direction which also agrees with the 
experimental data»^ A more detailed discussion of the 
momentum variation of the self-energy will be given in 
the final subsection. 

Because of the finite temperature in the ED/CDMFT 
calculation, it is not possible to identify spectral fea- 
tures at frequencies below the first Matsubara point 
(lj = 0.0314 for T = 0.01). Nevertheless, the doping 
variation of the pseudogap shown in Fig. [13] is found to 
be robust. In particular, it is clear that the pseudogap 
is directly linked to the resonance in — ImSjf(w) which, 
in turn, reflects the non-Fermi-liquid properties of the 
system. Since for 5 > S c ordinary Fermi-liquid behav- 
ior is established, it is evident that the pseudogap then 
vanishes. 

The above scenario is consistent with the fact that for 
a hole-doped Mott insulator the addition of electrons 
pushes the system closer to the insulating phase. This 
implies that spectral weight just above Ep must be re- 
moved and shifted towards the upper and lower Hubbard 
bands. This is precisely the effect induced via the large 
damping associated with the low-frequency resonance in 
-ImSj(w) and the positive slope of ReSx(w). 

According to this picture, the creation of holes in an 
electron-doped Mott insulator also moves the system 
closer to the insulating phase. Thus, spectral weight 
from states just below Ep must be shifted to the Hub- 
bard bands. As discussed in the next subsection, the 
ED/CDMFT results confirm this prediction. The X 
component of the self-energy along the Matsubara axis 
again exhibits non-Fcrmi-liquid behavior at sufficiently 
low electron doping. The extrapolation to real lo, how- 
ever, now reveals a resonance slightly below Ep, rather 
than above Ep as for hole doping. 



D. Electron Doping 

For completeness we discuss in this subsection the case 
of electron doping which differs from hole doping because 
of the second- neighbor hopping term t' . As a result of 
this interaction, the density of states shown in Fig. 1 is 
asymmetric, so that electron doping shifts the van Hove 
singularity away from Ef rather than towards it. Thus, 
the density of states is reduced and less steep. Fig. [TJ] 
shows the occupancies of the cluster molecular orbitals in 
the vicinity of the Mott transition induced via electron 
doping. Both these occupancies as well as the double 
occupancy shown in the lower panel exhibit hysteresis 
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FIG. 14: (Color online) Upper panel: Occupancies of cluster 
molecular orbitals (per spin) as functions of chemical potential 
IM, for U = 2.5, T = 0.01. For electron doping, the Mott 
transition occurs at about fi ~ 1.8 . . . 1.875. In the insulating 
phase 0.7 < fi < 1.8, nx = 0.5, nr = 0.75, and Um = 0.25 
(see also Fig. [3} . The arrows denote the hysteresis behavior 
for increasing vs. decreasing \i. Lower panel: Average double 
occupancy per site as a function of fj,. 



FIG. 15: (Color online) Upper panel: Imaginary part of 
Ex(iw„) as a function Matsubara frequency for electron dop- 
ing; U = 2.5, t' = -0M, T = 0.01. The long-dashed green 
curve for 8 = 0.12 approximately marks the transition from 
Fermi-liquid to non- Fermi- liquid behavior. The correspond- 
ing self-energy for hole doping is shown in Fig. [7] Lower 
panel: Comparison of damping 7x = —ImT,x(iu)„ —> 0) for 
hole doping (solid red circles) and electron doping (empty blue 
circles) . 



behavior for increasing vs. decreasing chemical potential, 
indicating that this transition is first order. Thus, this 
transition is similar to the doping-induced metal insula- 
tor transitions found within local DMFT for single-band 
and multi-band stystems i 80 ' 81 i 82 

Because of the lower and less steep density of states 
for electron doping, the low-frequency variation of the 
self-energy differs greatly from the hole doping case, as 
illustrated in Fig. [T5j Although there is again a clear dis- 
tinction between Fermi-liquid and non-Fermi-liquid be- 
havior, the transition now occurs at considerably smaller 
5. While for hole doping 5 C « 0.18 . . .0.20, for electron 
doping we find 5 C ~ 0.12. Thus, the Fermi-liquid phase 
is stabilized. 

To identify the pseudogap for electron doping, we eval- 
uate the cluster self-energy components via extrapola- 
tion to real frequencies. The upper panel of Fig. 1161 
shows at small to. In this case the non- Fermi- 

liquid properties give rise to a resonance in — ImSx(w) 
centered slightly below the Fermi level, indicating that 
the creation of hole states in an electron-doped Mott in- 
sulator implies a transfer of spectral weight from states 
near Ep to the Hubbard bands. Thereby the system is 
brought closer to the insulating phase. Accordingly, the 
real part of Ejf(w) exhibits a positive slope close to Ep. 



Its intersections with to + fi — can be used to define the 
pseudogap. For 5 = 0.08 the gap is found to be A sa 0.03, 
i.e., only about half as large as for the hole doping case 
shown in Fig. [T21 The lower panel of Fig. [T|5] shows the 
quasiparticle distributions obtained via extrapolation of 
the lattice Green's function components, Eq. (fl"T|) , to real 
to. The dominant feature at small to is the pseudogap in 
the X component, which is consistent with the behavior 
of displayed in the upper panel. 

The main difference with respect to hole doping, apart 
from the smaller size of the pseudogap, is the fact that 
this gap can be identified only in a very narrow dop- 
ing range. At electron doping larger than 0.08, the 
non-Fermi-liquid behavior is quickly replaced by ordi- 
nary Fermi-liquid properties. At smaller doping, spectral 
weight just below Ef is rapidly transferred to the Hub- 
bard bands, so that the pseudogap is superceded by the 
opening on the Mott gap. 



E. Phase Diagram 

In Figure[H]we have shown that the onset of non-Fermi- 
liquid behavior is shifted to smaller hole doping when 
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FIG. 16: (Color online) Upper panel: real and imaginary 
parts of Ejf(w) for electron doping obtained via extrapola- 
tion to real ui. The intersections of Re Ejc (a>) — [i with u> — et 
(straight line) provide the pseudogap. Lower panel: spec- 
tral distributions derived via extrapolation of lattice Green's 
function components G m (iiOn) to real u>. The average density 
corresponds to A{lu) = — ^Im [Gr(w) + Gm(lj) + 2Gx(w)]. 

U = 2.5 is reduced to U = 1.5 and when t' = — 0.3t is re- 
placed by t' — 0. Fig. fT5l illustrates the reduction of 8 C for 
(7 = 2.5 when hole doping is replaced by electron doping. 
A similar reduction is found for U — 1.5 (not shown). In 
Fig. [17] we collect these data and display the phase dia- 
gram of the present Hubbard model for electron and hole 
doping. At finite temperature the values of 8 C can only 
be determined within an accuracy of about ±0.02. For 
clarity, these margins are not plotted in Fig. [17] Despite 
this uncertainty, the results demonstrate several trends: 
for hole doping 8 C diminishes with decreasing U and when 
t' = — 0.3i is replaced by t' = 0. Moreover, for t' = — 0.3< 
the critical doping decreases when hole doping is replaced 
by electron doping. As pointed out above, the varia- 
tion of S c is surprisingly small, despite the rather large 
changes in U and t'. 

F. Momentum Variation 

According to the results shown in Fig.[S]the non- Fermi- 
liquid properties of the two-dimensional Hubbard model 
at low hole doping are mainly associated with the X com- 
ponent of the self-energy. Only very close to the Mott 
transition the M component begins to dominate since 
its imaginary part changes from ~ uj n to ~ l/w„. The 
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FIG. 17: (Color online) Phase diagram of two-dimensional 
Hubbard model calculated within ED/CDMFT for 2 x 2 clus- 
ters. The Fermi-liquid phase at large hole or electron doping 
is turned into a non-Fermi-liquid phase at small doping. The 
symbols for U = 1.5 and U — 2.5 indicate the approximate 
critical doping <5 C for t' = —0.075 (solid red circles) and for 
t' = (empty blue circles). The vertical line at 5 = marks 
the Mott phase at half-filling. The critical U indicated by x 
is about 1.4 for t' = and t' = — 0.3t. The long-dashed green 
line denotes the approximate lower bound of the non-Fermi- 
liquid domain. 

cluster components of the self-energy may be used to con- 
struct an approximate momentum dependent lattice self- 
energy by using the same periodization as in Eq. (|19[) for 
the Green's function. Thus^ 

E(k, u) = cv r (k)£rM + «Af(k)E M (w) + a x (k)E x (w) 

(20) 

where 

«r(k) = [1 + cosfc^Jfl + cosfcj,]/4 
a M (k) = [1 - cosfc K ][l - cosfc y ]/4 (21) 
ax(k) = [1 — cos/e a; cosfcj / ]/2. 

The k-resolved spectral distributions are then given by 
A(k,w) = --Im(w + /x-e(k) -E(k,w))" 1 . (22) 

An alternative is to periodize instead the cumulant 
matris^ M(ui) = l/[w + [I — E(w)] which can be di- 
agonalized in the same manner as the self-energy. Thus, 
the molecular orbital components of M(u>) are given by 
M m (uj) = 1 / [u> + [1 — E m (w)] and the momentum depen- 
dent lattice cumulant M(k, u>) can be derived from an 
expression analogous to Eq. ([20]) 

M(k,w) = ar(k)Mr(w)+aM(k)M M (w)+ax(k)J14(w). 

(23) 

The lattice self-energy in this approximation takes the 
form 

E(k,w) =u + /j,- 1/M(k,w). (24) 

In the upper panel of Fig. [15] we compare these two 
versions of E(k, ui) at u> = for 8 = 0.06 hole doping. 
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FIG. 18: (Color online) Upper panel: momentum variation of 
imaginary part of -Ex(u,k) at Fermi level. Solid red curve: 
expression Eq. (|20|l , corresponding to periodization of cluster 
self-energy; dashed blue curve: expression Eq. (|24|l . corre- 
sponding to periodization of cumulant; U = 2.5, 8 = 0.06, 
T — 0.01. Lower panel: azimuthal variation of — ImSx(w = 
0, k) for k = 7r(l — rsin<I>, 1 — rcos$) with r = 0.7, 0.8, 0.9, ap- 
proximately representing electron doping (e), half- filling (0), 
and hole doping (h), respectively. $ = 0,90 corresponds to 
k along XM and $ = 45 to the nodal direction FM. Solid 
red curves: periodization via Eq. (|20|l ; dashed blue curves: 
periodization via Eq. (|24|l . 



The real-w components T, m (tu — 0) are obtained via ex- 
trapolation from the first few Matsubara frequencies. At 
high-symmetry points both versions of ImS(k,w = 0) 
coincide. At general k-points, however, the cumulant ex- 
pression yields enhanced damping, in particular, between 
M and A, and along TM. The enhancement near A 
leads to an effective flattening of Im £(k, ui), which is also 
seen in the dual Fermion approach^ On the other hand, 
it is not clear whether this enhancement is partly an ar- 
tifact of the cumulant approximation since the damping 
at some points between A and M is even larger than at 
A. Also, damping near k w 2/3(jr, ir) in the cumulant 
version is almost as large as at A. At the present doping 
(5 = 0.06), the periodization of the self-energy according 
to Eq. (|20p is in better agreement with the dual Fermion 
approach (see Fig. 15 of RefiSi). 

The lower panel of Fig. [TH] shows the variation of 
— ImS(tj = 0,k) along k = 7t(l — rsm$, 1 — rcos<£>) 
where r = 0.7, 0.8, 0.9 is chosen to approximately repre- 
sent the Fermi arcs for electron doping, half-filling, and 
hole doping, respectively. Both periodization versions 
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FIG. 19: (Color online) Lattice spectral distributions A(k,uj) 
for three chemical potentials close to critical doping (fi = 0.3, 
upper panel), incipient pseudogap formation (p, = 0.4, middle 
panel), and complete pseudogap phase (fj, = 0.6, lower panel). 
Solid red curves (h): k x = k y = 0.36-7T near Fermi surface for 
hole doping; dashed blue curves (av): k x = k y — 0.5tt at 
center of Brillouin zone, corresponding to average density; 
long-dashed green curves (X): k = (tt, 0) = X; U — 2.5, 
T = 0.01. 



yield consistently larger damping along XM than along 
the nodal direction TM . The cumulant version implies 
overall larger damping, and, more importantly, less pro- 
nounced difference between TM and XM. Because of 
the substantial imaginary part of the self-energy at low 
frequencies, the Fermi surface in the present 2x2 clus- 
ter approach exhibits arcs rather than hole pockets. 40 
We note, however, that greater momentum differentia- 
tion obtained for larger clusters might lead to more pro- 
nounced anisotropy between the nodal and anti-nodal di- 
rections. In particular, this could yield smaller values of 
-ImS(w = 0,k) along TM (<f> = 45) than indicated in 
Fig. [H 
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FIG. 20: (Color online) Upper three panels: Spectral function 
A(\c, uS) along cuts 1, 2, 3 (from left to right); top panel: 
doping 8 = 0.17; middle panel: 8 = 0.14; U = 2.5, T = 0.005. 
Lower panel: Cuts through Brillouin zone corresponding to 
ARPES data in RefA Solid red curve: approximate non- 
interacting Fermi surface for hole doping. 



According to the collective mode in Im (oj = 0) (see 
Fig. [F3| , the anisotropy between the YX and YM direc- 
tions is even larger above Ep than for the u> — > limit 
shown in Fig. [TSJ Thus, the collective mode gives rise to 
a momentum and doping dependent particle-hole asym- 
metry. To illustrate this point, we show in Fig. [T9l the 
low- frequency part of the spectral distribution A(k, ui), 
derived via extrapolation of the lattice Green's function, 
Eq. (|19[) . at three representative points in the Brillouin 
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FIG. 21: (Color online) Same as Fig. [20] except <5 
(upper panel) and 8 = 0.08 (lower panel). 
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zone. Three doping regions can be distinguished: At 
/i = 0.3 close to optimal doping (5 = 0.18, upper panel), 
there is weak anisotropy since the system is a Fermi liq- 
uid throughout k space. Below critical doping (// = 0.4, 
5 = 0.12, middle panel), the spectrum in the anti-nodal 
direction at X shows clear signs of pseudogap behav- 
ior, while the one at k = 0.36(ot, 7t), i.e., near the nodal 
point of the Fermi surface for hole hoping, is still dom- 
inated by Fermi-liquid properties. At this k point, the 
coefficients in the momentum expansion, Eq. (|2U|) . are 
(ar, um , olx) — (0.50,0.09,0.41), indicating the rather 
large Fermi- liquid-like Y component. At the zone cen- 
ter these coefficients are (1/4, 1/4, 1/2). Finally, at even 
lower doping (/i = 0.6, 5 = 0.03, bottom panel), close to 
the Mott transition, the non-Fcrmi-liquid properties have 
spread across the entire Fermi surface, so that the pseu- 
dogap is observable along the nodal as well as antinodal 
directions. These results demonstrate the non-uniform, 
momentum dependent opening of the pseudogap as a 
function of doping. (Note that this behavior differs from 
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the opening of the Mott gap shown in Fig. [4] which in the 
present 2x2 cluster DMFT takes place simultaneously 
in all cluster components.) 

To analyze the particle-hole asymmetry observed in 
the recent ARPES data on Bi 2 Sr2CaCu208+«5 by Yang et 
alQ we have calculated the spectral distributions j4(k, u>) 
defined in Eq. (|22p , where the self-energy is obtained from 
Eq. (|20[) . The frequency variation of the X component is 
shown in Fig. [13] For direct comparison with the data we 
plot A(k, u>) along three cuts, as indicated in the lower 
panel of Fig. [20] Cut 1 corresponds to the nodal direction 
and has the lowest relative weight from E^(w), while in 
cut 3 the X component dominates. 

For large doping (Fig. [20] top panel), the system is a 
Fermi liquid. Thus, the spectral weight at all three cuts 
is largest at Ep and decays symmetrically for increasing 
and decreasing to. Below critical doping (middle panel), 
this particle-hole symmetry initially persists along the 
nodal direction, but gets weaker along cut 3. At <5 = 0.08 
(Fig. [2T] upper panel) , this asymmetry begins to extend 
to the nodal direction, until at S = 0.08 (lower panel) 
the particle-hole asymmetry is complete throughout the 
Brillouin zone. These spectral distributions reveal that 
the particle-hole asymmetry is a direct consequence of 
the pseudogap which gradually develops with doping at 
about 0.02 . . . 0.05 above Ep , and which is driven by the 
(■7T, 0) component of the self-energy. 

The momentum dependent opening of the pseudogap 
above Ep, and the particle-hole asymmetry caused by 
the collective mode seen in Im T,x (w) (see Fig. [TBI , are 
in excellent agreement with the ARPES data^! Although 
DMFT calculations for even larger clusters provide even 
better momentum differentiation, the present results for 
2x2 clusters reveal that spatial degrees of freedom give 
rise to dramatic new phenomena absent in a local descrip- 
tion, in particular, the resonance in the (ir, 0) component 
of the self-energy at small positive frequencies. It would 
be very interesting to check whether the dispersion of the 
position of this resonance with doping can be verified ex- 
perimentally. 

IV. SUMMARY 

The effect of short-range correlations in the two- 
dimensional Hubbard model is studied within finite- 



temperature ED combined with DMFT for 2 x 2 clus- 
ters. A mixed basis consisting of cluster sites and bath 
molecular orbitals is shown to provide an efficient and ac- 
curate projection of the lattice Green's function onto the 
cluster. The onset of non-Fermi-liquid behavior with de- 
creasing hole doping is evaluated for various Coulomb en- 
ergies, temperatures, and next-nearest neighbor hopping 
interactions. The self-energy component Sjjf = ( ff _o) ( w ) is 
shown to exhibit a collective mode above Ep which be- 
comes more intense close to the Mott transition. This 
resonance implies the removal of spectral weight from 
electron states above to Ep and the opening of a pseudo- 
gap. With decreasing doping the pseudogap opens first 
along the antinodal direction and then spreads across 
the entire Fermi surface. For electron doping, the res- 
onance of SxH and the corresponding pseudogap are 
located below Ep, as expected for the removal of hole 
states close to Ep. In the low doping range the density 
of states at the Fermi level becomes very asymmetric. 
Near the onset of non-Fermi-liquid behavior, Ep is at 
a maximum of the density of states. At smaller doping 
Ep moves into the pseudogap. This behavior leads to a 
pronounced particle-hole asymmetry in the spectral dis- 
tribution at intermediate hole doping, in agreement with 
recent ARPES measurements. The phase diagram shows 
that for hole doping 5 C « 0.15... 0.20 for various sys- 
tem parameters, i.e., near the optimal doping observed 
in many high-T c cuprates. The critical electron doping 
which marks the onset of non-Fermi-liquid behavior is 
systematically smaller than for hole doping. The Mott 
transition induced via electron doping exhibits first-order 
hysteresis characteristics. In contrast, within the present 
cluster ED /DMFT the hole doping transition appears to 
be continuous or weakly first-order at very low tempera- 
tures. 
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